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Abstract. We present observational evidence for the generation of Langmuir 
envelope solitons in the source regions of solar type III radio bursts. The solitons 
appear to be formed by electron beams which excite either the modulational instability 
or oscillating two-stream instability (OTSI). Millisecond data from the Ulysses Unified 
Radio and Plasma Wave Experiment (URAP) show that Langmuir waves associated 
with type III bursts occur as broad intense peaks with time scales ranging from 15 
to 90 milliseconds (6 — 27 km). These broad field structures have the properties 
expected of Langmuir envelope solitons, viz.: the normalized peak energy densities, 
^L/n e T e ~ 10 5 , are well above the modulational instability threshold; the spatial 
scales, L, which range from 1-5 Langmuir wavelengths, show a high degree of inverse 
correlation with ( Wz,/n e 7’ (f ) 1/ ' 2 ; and the observed widths of these broad peaks agree 
well with the predicted widths of envelope solitons. We show that the orientation of 
the Langmuir field structures is random with respect to the ambient magnetic field, 
indicating that they are probably isotropic structures that have evolved from initially 
pancake-like solitons. These observations suggest that strong turbulence processes, such 
as the modulational instability or the OTSI. stabilize the electron beams that produce 
type III bursts. 



1. Introduction 


The production of solar type [If radio emission involves the generation of high 
levels of Langmuir waves and their subsequent conversion into electromagnetic radiation 
at fpe and 2f pe (f pf = [n e e 2 /irmg) { l 2 is the electron plasma frequency, where e, n e , and 
nig are the electron charge, number density, and mass, respectively). Electron beams 
propagating outward from the sun amplify Langmuir waves through the bump-on-tail 
instability [Bohm and Gross , 1949J at the resonant frequency u>l = kiv b , where k £ 
is the wave number of the Langmuir wave, v b is the speed of the electron beam, 

= •*£.+ 3A: 2 vf e , u pe = f P g, and v Te is the electron thermal speed. Early theoretical 
models describing the injection of electron beams into the solar atmosphere predicted 
that excitation of Langmuir waves would extract all the streaming energy from the 
electrons within 100 km or less f Sturrock , 1964: Tsytouich, 1970], whereas in situ 
detection of type III burst associated Langmuir waves and electron beams [Lin et al., 
19S6; Kellogg et al., 1992; Reiner, Fainberg and Stone, 1992; Thejappa et al., 1993a, b; 


Thejappa. Wentzel and Stone. 1995; Thejappa. Stone and Goldstein. 1996: Thejappa and 
Mac Do wall. 199b] indicate that the electron beams reach 1 Astronomical L nit (AL) or 
more. Consequently, the beam cannot be in continuous resonance with Langmuir waves 
while propagating from the corona to and beyond 1 AL . Some process must remove 
Langmuir waves from resonance with the beam by, for example, reducing k^ through 
uvak turbulence processes, such as induced scattering or electrostatic decay [see, Kaplan 
and l sij/onch. 19i3] or by increasing k/ through st rong turbulence processes such as 
I lie modulatioual instability or oscillating-two-st ream instability (OTSI) [Fapadnponlos, 
Goldstein and Smith. 1971; Smith. Golds! t in and Tnpadopoulos, 1979: Goldstein. Smith, 
and Tnpadopoulos, 1979], 

Beam stabilization occurs via weak turbulence processes, when the rate of removal 


of Langmuir waves from resonance with the beam exceeds tin* growth rate of the 
instability. In contrast, the wavepackets generated by strong turbulence exert a 
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pond< lomoti v< fou, on the plasma, making it spatially inhomogeneous. 1,1ms causing 
a loss of coherence and decoupling the wavepacket and beam. In the first application 
of strong turbulence theory to the type III burst problem, Papadopoulos, Goldstein 
and Smith [1974] used the term oscillating two stream instability, or OTSI, to describe 
the initial strong turbulence processes that form the small-scale soliton-like structures. 
During later stages of evolution, these solitons collapse [ Zakharov 1972: Nicholson et al. 
19(3; Kellogg et al. 1992; Shapiro and Shevchenko , 1984; Robinson 1997], 

The linear regime of the OTSI or modulational instability is characterized by the 
formation of an envelope soli ton which is the envelope of the high frequency Langmuir 
waves trapped inside a self-generated density cavity. The envelope soliton is also 
referred to as an oscillating soliton [ Kuznetsov . Rubenchik and Zakharov , 1986], because 
it contains oscillations of definite frequency and wavelength. The size of the envelope 
soliton, which ranges from 1 to several Langmuir wavelengths is determined by the peak 
intensity of the trapped waves: the higher the field strength, the narrower the width. 

In previous studies, we evaluated the emission mechanisms at f pe [ Thejappa et 
al.. 199-3a.b] and at 2/ p( . f Thejappa . Stone and Goldstein , 1996] using in situ wave 
data associated with several type III bursts from the Unified Radio and Plasma 
Wave Experiment (LRAP) on Ulysses. Our main emphasis in those studies was to 
evaluate various type III emission mechanisms. In the present paper, we reanalyze 
the high resolution observations of Langmuir waves that occur as the broad intense 
field structures with durations ranging from I") to ')() milliseconds. We argue that 
these field structures are probably quasi-stable Langmuir envelop*' solitons generated 
by the modulational instability. Our conclusion is based on the fact that the peak 
intensities are well above the modulational instability thresholds, the spatial scales 
(1-5 Langmuir wavelengths, X L ) are consistent with the widths computed for envelope 
solitons, and, most importantly, the widths of these envelopes vary inversely with their 
iiil.rnsiUrs. 
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Tor most of the* Langmuir wave broad peaks, the normalized energv densities 
^ ) are slightly less than or equal to the corresponding values of . where 

fL- — tB/(m r _c) is the electron cyclotron frequency. A careful examination of the data 
has revealed that the ambient magnetic field B is oriented randomly with respect to 
the antenna direction, which is the probable direction of the high resolution events of 
Langmuir wave field structures. This implies that even though magnetized Langmuir 
waves initially form pancake-like structures with transverse spatial scales 5 X much 
larger than the longitudinal ones (S), i.e., S± S , they evolve into more isotropic 

structures with S± ~ 5. This isotropization occurs because the perpendicular spatial 
scales decrease more rapidly than do the parallel spatial scales so long as the inequality 
holds [lirasnoselskikh and Sotnikou , 1977], YVe also show observations 
of weak spectral enhancements of ion acoustic waves at about ~ 100 Hz in association 
with the envelope solitons. YVe discuss the implications of these observations for the 
stabilization of the electron beams that excite type III radio bursts. 

In section 2. we present the observations. Section 3 contains a solution of Zakharov’s 
equation (unmagnetized case) in the form of an envelope soliton and compares it with 
observations. The magnetized Langmuir waves are shown to form planar solitons that 
evolve into isotropic solitons and eventually collapse. In section 4. we discuss the 
implications of these observations for type III burst beam stabilization and the possible 
coexistence of electrostatic decay with modulationa! instabilities. In section 5. we 
present the conclusions. 

2. Observations 

VV'e concentrate primarily on three local type 111 events and their associated in 
situ waves, observed by the UHAP experiment on II December 1990, 22 February 
1991 ( identified by Reiner, Fainbenj and Stone, [1992)), and 7 March 1991 (identified 
by The jap /mi rt at., (1993a)). In addition to a variety of wave data from the Radio 
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Astronomy Receiver (RAH). th«' Plasma Frequency Receiver (PFRj. the Wave Form 
Analyze! (WPA), and the Fast Envelope Sampler (FES) of the FRAP experiment [Stone 
tt «/., 1992], we use magnetic field ( B ) data from the fluxgate magnetometer [Balogh et 
al., 1992], and the electron density (n t ), electron temperature ( T e ). ion temperature (T,), 
and solar wind speed (V' JTU ) data from the Solar Wind Plasma Experiment SWOOPS 
[ Bame et ai, 1992]. These data were provided by the National Space Science Data 
Center (NSSDC). 

In columns 1 , 2 and 3 of Figure 1 , we present the peak electric field signals 
corresponding to local type III events of December 11, February 22 and March 7, 
respectively. The first and second rows show the time profiles of the type III bursts 
at ~ 2 f pe and ~ / pe , respectively, whereas the bottom row shows the corresponding 
Langmuir wave electric field signals. It is clear from this figure that the type III burst 
time profiles are generally characterized by a steep rise followed by a smooth and 
slow decay. Since the second harmonic emissions do not have exactly the same time 
development as that of fundamental emissions, it is very likely that the high frequency 
electromagnetic emission contains the first as well as the second harmonics. It is clear 
from the bottom row that the Langmuir wave electric fields are very bursty in nature. 

1 heir peak amplitudes are 3 to 4 orders of magnitude larger than the respective type 
III burst electric field amplitudes. The type III burst electric field measurements are 
from the RAR with a time resolution of ~ 128s and a high sensitivity of ~ 30nVHz -l/ ^ 
al 10 kHz. The Langmuir wave data are from the PFR with a time resolution of ~ His 
(maximum value in 16 s windows) and a rms sensitivity of ~ IOOiATIz -1 ^ 2 at. 10 kHz. It 
is clear from this figure that the intense Langmuir wave' activity coincides with the peak 
of the type HI burst emission. In Table 1, we list all relevant solar wind parameters, 
such as the electron density, n K \ the electron temperature, the ion temperature, T,\ 
the solar wind speed, and the magnetic field, B. 

figures 2a, I), and c show the high time resolution snapshots of Langmuir wave 



,rrlr "‘ t * M v**lop**s associated with the type III bursts of December 
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Tebruary 22, 1991, and March 7. 1991, respectively. These data are obtained by the 
FES, which is capable of resolving the field structures with time scales as small as 
one millisecond (for a detailed description of the instrument, [see Kellogg el al., 1992; 
Theissen and Kellogg, 1993]. As seen in Figure 2, the Langmuir wave events show 
tandom fluctuations in the beginning and a broad intense peak in the middle, after 
~ 500 milliseconds. Step-like structures are often present during the rising phase of 
these broad peaks. The asymmetric shapes and enhanced noise levels during the decay 
phase of these events are due to the 32-dB attenuator, activated by the intense central 
peak (but cf. panel b). The role of the attenuator is to extend the dynamic range of the 
system, i.e., whenever the signal rises to saturation, the attenuator is switched into the 
signal stream, which remains connected for the remainder of the event cycle. 

We have examined the orientation of the solar wind magnetic fields with respect 
to the X-antenna. During most of these FES events, the magnetic field is oriented 
randomly with respect to the X-antenna. Because the FES captures the most intense 
signals during each 30 minute interval, the observed events probably the electric field 
structures aligned along the antenna direction. This suggests that the Langmuir 
wave electric field signals of this study are also oriented randomly with respect to the 


magnetic field which implies that the Langmuir envelopes are approximately isotropic. 
In Table 2. we present the summary of these high time resolution observations. These 
include: (I) the peak electric field intensities of the broad central field .structures. K L . 
(_) the normalized Langmuir wave peak energy densities, I n,.T r ) = c„ /'.'//( 2n, /',) 
(<„ is the dielectric permittivity of the free space). (3) the measured 0.2-power duration 
of broad peak, (t 0 , 2 ) in units of milliseconds, (d) the 0.2-power widths of the broad 
peaks, L 0 i = t 0 ,iV, w , and (.5) the predicted 0.2-power widths of envelope solitons S 0 . 2 
calculated using the observed and X t) (see, section. 3). 
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l.lic square root, of the corresponding inverse normalized energy densities. [{n r [',.)/ 
is presented. As seen from this figure, £ 0.2 increases with implying 

that the more intense the peak, the narrower the width. We compute the correlation 
coefficient in this case as ~ 0.9. Such an excellent correlation indicates that these 
observed field structures most probably correspond to Langmuir wave fields trapped 
inside self-generated density cavities. 

In Figures 4a, b and c, we plot the low frequency electric field spectra observed 
by the WFA in the frequency range 0 — 448 Hz, as well as high frequency electric field 
spectra observed by the PFR (0.57 — 35 kHz), during I- or 2-minute intervals containing 
the times of the FES events. The dotted lines correspond to the instrumental threshold 
of WFA. The prominent spectral peaks at ~ 10 kHz correspond to Langmuir waves, and 
weak spectral enhancements at ~ 100 Hz correspond to ion acoustic waves. 

In the following sections, we examine whether the prominent broad central 
peaks are quasi-stable Langmuir envelope solitons and discuss their possible 
generation mechanisms. We discuss briefly whether the observed electric field spectral 
enhancements at low frequencies during the intense Langmuir wave envelopes are due to 
the coexistence of the electrostatic decay (weak turbulence) and modulational instability 
(strong turbulence) processes in some of the type III burst source regions. 

3. Modulational Instability and Langmuir Envelope Solitons 

Neglecting the magnetic field, the dispersion relation ot Langmuir waves can be 
written as: 

U '*Z = » + ; ^7. (1) 

where /* 7V is the electron thermal speed. Whether the beam excited Langmuir wave 
packets broaden due to dispersion (weak turbulence processes), or, shrink due to 
ponderoniotive effects, is determined by a threshold (neglecting dissipation) [Sa<]dn t\ 
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1979; Shapiro and Shevchenko, 198-lj: 


~^3(Ak L A D )\ (•; 

n e T e '■ 

where Ak^ is the average spread in wave numbers of the Langmuir wave packet. This 
is ie!ated to the typical values of the velocity dispersion in the beam, Av h , and the 
observed number of linear growth times. A\ of Langmuir waves before onset of the 
modulational instability as [cf. Lin et al. 1986]: 


Aki _ Av b In 2 

~kT ~ ~^2N' 


A, the number of linear growth times can be estimated from N ~ Inp, where p is ratio of 
the peak electric field amplitude of the broad peak (Figure 2) to the thermal background 
which is approximately an order of magnitude lower than the instrumental background. 
Here, vve note that no self-consistent procedure is available for incorporating the 
collisionless dissipation in the threshold criterion (2), however, it is most likely that 
the dissipation processes slightly increase the modulational instability thresholds in 
comparison to those estimates obtained from Eqn (2). In the present case, the observed 
ratio of p ~ 10 3 corresponds to .V ~ 7. So, for Av h /v h ~ 0.15 [Lin et al. 1986] and 
•V ~ "• vve obtain Ak L /k L ~ 7 x 1(T 3 . For typical beam speeds, v b ~ 6 x 10 7 ms' 1 


[Erg an et al.. 1998], and the observed values of the electron plasma frequencies f pe 
(Table I), we obtain k L ~ u>/e 6 - 1.2 x HT* m" 1 . - [.[ x 10" 3 m" 1 and ~ U x 10" ! 


, - 1 


for the December I I. February 22, and March 7 type HI events, respectively. Tin 
corresponding wavelengths are \ L - >ir/k,, ~ 5. 5.7 and 5.7 km. ( sing the estimated 
values of Ak L /k L ~ 7 x l() -3 , we find Ak { , ~ 8.1 x I0 _,i m -1 . — 7.7 x |()-' ! 

~ x 10 ,J m -1 for the three events. The observed values of Debye lengths (\ 0 ) 
corresponding to these events are 17, 17 anti LI meters (see. Table I). By using these 
values of Ak ts and X D , the modulational instability threshold W tl J(n r T i: ) is estimated 
as ~ 6 x 10-*, 5 X 10-* and 5 x 10~« for the three events. Here, vve note that the initial 
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l><-am excited Langmuir wave numbers are k L X a < am | the corresponding 

group speeds are l g = :f(k L X d )v Tk - (0.04 - 0.06)1*, - c,. where m ( and m, are the 

electron and ion masses, and V Tt and c, are the electron thermal speed and ion acoustic 
speed, respectively. 

From Table 2, it is clear that the observed normalized peak energy densities 
^L/(n e T e ) significantly exceed the modulational instability thresholds. This implies 
that the ponderomotive force effects that drive the modulational instability are stronger 
than dispersive effects. The ponderomotive force arises from spatial gradients in wave 
intensity which affects both electrons and ions equally causing them to move toward 
the intensity minima, thus forming a ripple, 5n e , in the plasma density. The waves get 


trapped in those regions where the density is low because, from the Langmuir wave 
dispersion relation, waves with large k L exist only where uj pe is small. The wave trapping 
fuither enhances the intensity in the regions where it was already high, thus causing a 
ripple in the envelope to grow. For the modulational instability to occur, it is essential 
that the nonlinear frequency shift Su> L , responsible for the shrinkage of the wave packet, 
should have a sign opposite that of the group dispersion dv g jdk L which tends to spread 
the wave packet. Thus, after the modulational instability has grown somewhat, the 
Langmuir wave electric fields in real space consist of wave packets with group speeds less 
than ~ c 3 . These wave packets will begin evolving into a series of Langmuir envelope 
soli tons with group speeds less than - c, [Degtyarev. VV ikhankor and Rudakor. 1971; 


\n'holson ft al |<)7Sj. 

3. 1 Envelope Solitons in an Unmagnetized Plasma 


I’he set of equations describing the modulational instability processes, which 
UK lude envelope soliton formation are the well-known Zakharov equations [Zakluiroi 
I972J. which, in one dimension, take the form 
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r)E d 2 E 

‘aZ + W = nE W 

d 2 n d 2 n d 2 \E 2 \ 

dr 2 dz 2 dz 2 ^ 

Here we use the dimensionless variables defined by Nicholson [1983]: 



where 7 e and 7 ,• are the specific heat ratios of electrons and ions, respectively. The 
solutions of equations (4) and ( 5) can be found in terms of envelope solitons moving with 
constant speeds V < V g , where V 3 is the group velocity of the Langmuir waves. Several 
authors [see. Rudakov . 1973: hingsep. Rudakov and Sudan. 1973: Nishikawa. Hojo and 
Mina. 1974: Schamel. Yu and Shukla. 1977: Kuznetsov . Rubenchik and Zakharov , 1986] 
have obtained the envelope soliton solutions by solving equations ( 4 ) and (5) under 
different approximations. These authors either have neglected the ion temperatures, 
i.e., T, ~ 0. or have solved the nonlinear Schrodinger equation (obtained by neglecting 
the first term in equation (5)) for stationary solitons. In this paper, our motivation is to 
obtain the soliton solutions for realistic solar wind conditions and compare them with 
observations. The Mach number of the envelope soliton is defined as: 



where the ion sound speed c a is defined as: 


( 11 ) 


r ; = 


7 r l r + 7 , 7 ', 


m, 


( 12 ) 
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VI 1 1. hoi it loss of generality, one can take y K = 1 and 7, = 5/2. 
The solutions for equations (4) and (5) can be written as: 

E = E(£) exp (— i0(r, £)) 

n = n(0. 


( 13 ) 

(14) 


where <f = r - A/r and is the phase of the envelope. Because, these are the solutions 
of solitons moving with constant speeds V, the variables E and n depend only on 
(* - Mr) but not on the spatial coordinate, z and time, r, separately. Equation (4) can 
be separated into imaginary and real parts: 


and 


h ,dE 
-M- + 


[dr M dt\ 


'd 2 4> L ,dEd<t > 1 n 

E w +2 *dtr° 


d 2 E 

d? 


~E(p]=nE. 


^From equation (15), the phase of the envelope can be obtained as: 

= -ftr + 

where Q is some constant frequency. By using d 2 /dr 2 -* Md 2 /6£ 2 and d 2 /dz 2 
in equation (5). we obtain 


n = 


(15) 

(16) 

UT) 

d 2 /d ( 2 

( 18 ) 


i - \n ' 

By substituting n and d>(r,£) from equations ( IS) and (17) in equation (16). one obtains: 


8 2 E 

d? 


= -E 




E' 


i - \n 


14) 


liy integrating once, this equation becomes: 


M) 


M M 2 

{} + — ■ — -f- 


E‘ 


4 2( l - M 2 ) 


( 20 ) 


Writing E(0 = E 0 f(Z) and 


!! = -—- 


Hi 


I 2( I - M 2 ) ’ 


( 21 ) 



equation (20) becomes 



Ej 

2(1 - \P) 


*(£) a [i -f(O a J. 


(22) 


where Eq is the normalized peak intensity of the envelope, and e(£) is some auxiliary 
function. By making the change of variables 


£ = 


y2(l-A/’) r , 

Eo * 


(23) 


equation (22) becomes 

(^P) I=£(f ' , ’[ 1 - £( f ' )2 J- (24) 

The solution of this equation, which decreases as if -> ±oo can be written as [Tsytovich, 
1995]: 

c = cosh -1 <f\ (25) 

Changing back to the original variables (equation 23), equation E(£) = 

£oe(£)exp( — i0(r, £)) becomes 


E = Eq cosh 


-l 


Eot 


■ (MZ n \ 

exp — i I — J2r I. 


72(1 - W ) 

( sing equations (6) — (10), this soiiton solution can be expressed as 

/ j j/ ^ \ 

El = fi’iocpsh -1 ^ — — — j cos(u,'yt — k 0 z), 


(26) 


(27) 


uhete Ei is the amplitude of the envelope and Eiq is its peak amplitude. These 
envelope solitons can be stable for soiiton speeds l < (l /y/$)c a [see, Vladuntror rt ai. 
1995]. 

Let us discuss the significance of each term in expression (27) in the light of the 
present observations. For example, ujy is the nonlinear frequency shift, which is related 
to the observable parameters as: 
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I'or envelope solitons with spools [ - V a . tho nonlinear frequency shift, bec omes: 


W- 


W 


l*D 




Sn e r r 


(29) 


i.e., the nonlinear frequency shift a/.v is due to the dispersion (first term in (29)) as well 
as the ponderomotive effects (second term in (29)). These two effects oppose each other, 
i.e.. have opposite signs. These types of nonlinear frequency shifts are very difficult to 
measure because they are too small (for example, in the case of a stationary solitons 
with M -f 0, u;/v/u; pe ~ -W L /8n e T e , which is ~ 10~ 6 for the events of the present 
study). The central wave number of the envelope soliton k 0 (equation 27) is defined as 

M (r)m t \ */ 2 1 


k = — (T2i \ 1/2 1 
3 \ m, / X D ' 


(30) 


i.e., k 0 ~ k L for V - V 3 . The soliton width 5 is related to normalized peak energy as 
(equation (26)): 


Ad. 


(31) 


s- = 'h] l/2 

1 1 - A/ 2 W L J 

i.e., 5 oc A o (^L/n e r e ) 1/2 . This relationship between the soliton width and the 
peak energy density can be observational ly verified. For example, in Table 2. the 
observed values of W L /{n e T e ) and the 0.2-power widths L Q2 of the broad peaks are 
given. In Figure 3 we plot the measured 0.2-power width of the broad peaks L 0 . 2 versus 
( £./'** T e ) , which shows that higher the ( H'£,/n s 7’ e ) 1//2 , narrower the width £ 02 . 

We calculate the correlation coefficient between these observed quantities as ~ 0.9. 

I Inis, the observations agree very well with the expected scaling for Langmuir envelope' 
solitons (equation 31). This implies that the Langmuir wave* bursts captured by the 

L L.S instrument are most probably the fields trapped inside the self-generated density 
cavities. 

We can also compute the predicted widths of envelope solitons using equation 
(31). since all the quantities on the right hand .side (IMIS) are measured except the 
Mach number M, which can be taken as - 1/^3 corresponding to the stable Langmuir 
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envelope solitons. For tin? measured temperatures of electrons (T„) and ions (T t ) as 
given in Table 2, we obtain r, = ( lr T r + y.T.) /T, - 2.25, 1.33 and ~ 1.37. for December 
11. February 22. and March 7 events, respectively, where we have assumed ~/ e = l 
and 7, = 5/2. We denote the envelope soliton width, where the field is 0.2 E^q as 
S 0 . 2 - Therefore, from equations (26) and (27), we obtain cosh“ l (.5o. 2 /5’) = 0.2 for 

(‘S.i/'-’) = ~-29. i.e., for t = 0 and r ~ So. 2 - By using this, we can write the 0.2-power 
width of the envelope soliton as: 

c A ( V n e T'\ l /\ 

A °' < 32 > 

Thus, we obtain: S 0 . 2 ~ 6.7 (n e T e /W L ) 1/2 X D . ~ 5.1 {n e TJW L ) l/2 \ D , and ~ 

5.2 {n e T e /\\ l) \ d , for type III events of December, February and March, respectively, 
where we have used XI = V/c, - 1/^5. Using (WT/n c T e ) and from Tables 2 and 
1, we estimate the expected values of 0.2-power widths of envelope solitons 5o. 2 - The 
results are shown in Table 2 where the predicted value S 0 . 2 as well as the measured 
value Lq. 2 , are included. The 5o. 2 and L 0 . 2 agree remarkably well, except, perhaps, for 
the second FES event of December 11 and first FES event of March 7 type III bursts. 
From Fig. 3, one can see that these two points lie well below the fitted line. 

Using Lq .2 and the wavelengths of the Langmuir waves, we can also estimate 
the number of Langmuir waves trapped inside these envelope solitons. For example, 
the envelopes of December II contain 5. 3 and I Langmuir waves: the February 22 
envelopes contain 3 and 5 waves; and the March 7 envelopes contain 2. 1. 2 and 2 waves, 
respectively. In terms of Debye lengths. \ 0 . the widths of envelopes of December I I. 
February 22 and March 7 are equal to 1320, 700, 1230, 1160, 160. 710. 480. 830 and 730. 
respectively. Thus, we conclude that the observed broad field structures are Langmuir 
envelope solitons because the measured spatial scales are in excellent agreement with 
the predicted widths, .\ J2 and because of an excellent correlation between the L {) 2 and 
as expected for envelope solitons. 
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3.2 Langmuir Solitons in a Weakly Magnetized Plasma 

fn the case of type III bursts, although the ambient magnetic field B is very weak, 
it guides the electron beam from inner solar atmosphere to several AU in the solar 
wind, thus the beam velocity v b is almost parallel to the ambient magnetic field, B. The 
weak magnetic field also determines the evolution of the beam excited Langmuir wave 
spectrum, which is forward peaked with a growth rate of (Kaplan and Tsytovich, 1973): 
lb/wp e = n b /n e (v b j Av b ) 2 cos 9, where n b and n e are the electron densities of the beam 
and the ambient plasma, respectively, v b is the beam speed, Av b is the spread in the 
beam speed and 9 is the angle between xy, and k i. Since this growth peaks at 9 ~ 0, the 
initial spectrum of beam excited Langmuir waves is one-dimensional in nature. 

When the ambient magnetic field is taken into account, the dispersion relation of 
Langmuir waves becomes: 


= U 


P« 


, , 2-, n; ki 
1 + 2^°' + ■ 


(33) 


where is the component of the wave vector perpendicular to the ambient magnetic 
field. In the static approximation, the envelope equation for the amplitude of the 
magnetized Langmuir waves in dimensionless variables has the form [ Dysthe et al., 
197SJ: 


d 2 * , d 4 * t 3 d 




dz 


2 , 


d* \ n 

d =0 ' 


(34) 


1 dz 2 ' dz* 1 ^ * ■ dz 

where is the envelope of Langmuir waves, t(( = 1 - ^ and t L = 1 - j ■ Bv a 

muiH'iiial method, a pancake shaped solitou symmetric around an axis parallel to B j 

found for this equation by Pelviashvili (1973}. He has shown that the thickness of this 
soliton is: 

S = A 0 /( W L /n r T r ) {, \ (35) 

whereas the radius of the pancake is 


- s \ =S(iX:/~'l 


(36) 
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! Inis, m the early stage of soliton formation, as long as 1 >*/-£ » k[Xj } . the transverse 
dimensions of the pancake-like solitons (S L ) are much larger than the longitudinal ones 
(5), forming highly anisotropic, elongated, dipole field structures. Dysthe et al [1978] 
have found an analytical one-dimensional soliton (planar) for equation (34) as: 


E l = E l osech S~ l (z±sf + 2e ll 


1/2 S 


L 


2|ex| 


(37) 


where 5 is the thickness of the pancake as given by equation (35), the local magnetic 
field is assumed to directed along z-axis. The stability of these solitons against the 
transverse perturbations has been studied by several authors [Rowland, 1985; Kuznetsov, 
Rubenchik and Zakharov, 1986; Hadzievski et al., 1990; Hadzievski and Skoric, 1991], 

The linear growth rate of the soliton instability is [Kuznetsov. Rubenchik and Zakharov, 
1986; Hadzievski et al., 1990]: 


7/^j 


pe 




(38) 


where ^(j) ^ nU j s Rjemann s zeta function. This shows that a strong enough 
magnetic field B with 


9.1 W, 

> 0.43— 

w'pe n f T K 


(39) 


can stabilize the linear instability. However. Hadzievski et al [1990] have shown that 
the instability reappears if one takes into account the next term in the expansion in 
transverse wave numbers k x , and Langmuir collapse will ensue. 


As discussed by Krasnosdskikh amt Salnikov [ 1977 ]. even though for 9 .,. /„ 


*P* 


>> 


k,\n. the transverse dimension of the planar soliton ,s' x is larger than the longitudinal 
one S. and A'i will decrease rapidly in comparison with .V. until the solitons become 
isotropic. Uius, the Fast Envelope Sampler appears to capture the solitons at this 
isotropic stage of their evolution when k L X D is approaching Q r /u. w . and the solitons 
evolve independently of the local magnetic field. These solitons should collapse 
< ventually. I hose observations do not confirm the simulation results of Hadzievski rl al. 
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[1990], which showed that, tin* form of the cavities and trapped energy depend on the 
magnetic field strength. 

If the soliton conserved energy, Krasnoselskikh and Sotnikov [1977] showed that 
they could not construct a self-similar solution to collapse in a magnetic field. If the 
soliton radiated away the electrostatic energy, however, they could get such a solution. 
Rotvland [198-5], using computer simulations, showed that if Wi/(n t T e ) < (f> e /u/ pe ) 2 , 
planar solitons with a finite transverse scale would radiate the wave energy before 
collapsing. This agreed with the earlier theory of Krasnoselskikh and Sotnikov [1977]. 
This process is termed as weak collapse by Kuznetsov and Turitsyn [1990] and Hadzievski 
et al [1990]. When the most of the energy is absorbed by the electrons, the collapse 
process is called strong collapse [Kuznetsov, 1996]. 

4. Discussion 

We have analyzed the ~ 1 millisecond observations of the Langmuir waves which 
form broad intense peaks with spatial scales ranging from 480 to 1320 Xp. Our 
analysis indicates that the peak intensities of these envelopes exceed the threshold for 
modulational or OTSI instabilities. The widths of these broad Langmuir field structures 
correlate well with (n e T e / as expected for envelope solitons. This implies that 
the density cavities associated with these, solitons are most probably generated by the 
Langnmir waves themselves due to their ponderomotive force. The predicted widths of 
envelope solitons calculated using the observed intensities of the* broad Langmuir field 
structures agree remarkably well with the observed widths. These findings strongly 
suggest that the broad Langmuir wave peaks correspond to envelope solitons generated 
by the modulational or oscillating two stream instabilities. This implies that strong 
turbulence processes play a significant role in the stabilization of the electron beams 
which generals' the Langmuir waves. 

In earlier work [The jap pa et at., 1993a], we discussed the role of the electrostatic 
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<l<‘< ay instability which is a weak turbulence process in the development of type III 
bursts. During this process, the beam-excited Langmuir waves are expected to decay 
into daughter Langmuir waves and ion acoustic waves. As shown in Thejnppa et al., 
[199.3a], the threshold condition for electrostatic decay is satisfied for all broad peaks 
shown in Figure 2, implying that this process should co-exist with strong turbulence 
processes in the source regions of some of the type III events. One can estimate the 


maximum possible Doppler-shifted frequencies f s for the daughter ion acoustic waves 
by using the relation: /, ~ k s (V sw + c,)/(27r), where k s ~ 2 k L is the wave number of 
the ion acoustic waves. Using the observed values of V 3Wl c, and k L (see, Table 1), we 
obtain the predicted values for f s as 99.3 Hz, 69.2 Hz, and 86 Hz for the Dec. 11, Feb. 
22 and March 7 events, respectively. One can also estimate the expected ion acoustic 
electric field values, E„ by using the energy conservation arguments. In the present 
case. E s values lie in the range (1.4 - 3.8) x 10~ 5 Vm" 1 . (5.3 - 6.7) x 10" 6 Vm" 1 and 
(0.9 -2.7) x 1 0- 5 Vm" 1 for the events of Dec. 11, Feb. 22. and March 7. In our previous 
study, we failed to find any electric field enhancements which might correspond to these 
predicted ion acoustic wave signals. We therefore concluded that although threshold 
conditions were satisfied, the electrostatic decay process was not evident. 


In this study, rather than using the time profiles as in Thejappa et al. [1993a], 
we analyzed the spectral data of peak electric fields obtained by the VVFA during 1 
to 2 minute intervals covering the FES events. In Figures -la. b and c, we present the 
low frequency electric field spectra from U FA as well as tin- high frequency spectra 
bom the PER. It is clear from these plots that for most, of the cases, the low frequency 
electric field spectra show weak enhancements at - 100 11/ coinriilent with the passage 
of modulationally unstable soli tons. We interpret these simultaneous occurrences of low 
frequency electric field signals and Langmuir waves as evidence of the coexistence of 
weak and strong turbulence processes in the type III burst source region [Thejappa and 
Mar l)a wall, 1998]. This suggests that during some type III events the electrostatic decav 
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instability can ivmov, Langmuir wav,* from resonance with electron Loams. l i.is loads 
to the accumulation of Langmuir waves at long wavelengths, forming the so-called weak 
turbulence Langmuir condensate, which eventually becomes modulationally unstable. 

There are at least two reasons for not detecting low frequency wave enhancements in 
association with solitons for some of the events. First, the size of the solitons is < 40 km. 


well below the coherence length of v h / lL > 1000 km. Therefore, the resonances required 
by weak turbulence processes cannot occur. Second, the threshold for electrostatic decay 
is at least two orders of magnitude greater than that for modulational instability (see 
Sec. 4), thus electrostatic decay should not be excited, at least until the modulational 


instability has saturated. 


High frequency ion acoustic-like waves [Gurnett and Frank. 1 978] can play the role of 
small scale (k, >> k L ) density fluctuations: which can cause the nonresonant anomalous 
absorption of Langmuir waves due to scattering [Goldman and DuBois. 1982: Russell 
and Goldman. 1983]. This can also stabilize the electron beam. The observed lack of 
high frequency ion acoustic-like noise between 0.5 to 3 kHz during intervals of intense 
Langmuir wave activity was interpreted in terms of this process. As shown previously 
[ Thejappa et. al. 1993a]. the ion acoustic-like waves are not associated preferentially with 
Langnnm waves, i.e., Langmuir wave growth is not suppressed by scattering on density 
fluctuations associated with high frequency ion acoustic-like waves. This implies that 
anomalous scattering of Langmuir waves by small scale density fluctuations is probably 
not a significant stabilization mechanism for type HI burst exciters. 


5. Conclusions 


Lsing Tlvsses data from the 


TRAP 


experiment we have presented 


observational 


evidence for the Langmuir envelope solitons, generated by the modulational instability 


or OTS'l in the source regions of solar type HI 
observations show that the type III associated 


radio bursts. The high time resolution 
Langmuir waves occur as broad intense 
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[><viks With time scales ranging from io - '10 ms (equivalent to the spatial scales of 
f> “ 27 km) ' We have identified the broad peaks as Langmuir envelope solitons based on 
(1) M L /n e T e of the broad Langmuir peaks are well above the modulational instability 
threshold of ~ 10" 5 , (2) the spatial scales of these localized field structures, which 
range from 1 to 5 Langmuir wavelengths, show a high degree of inverse correlation with 
(Wl/I'T') 1 '*, as expected of envelope solitons. and (3) the observed widths of these 
broad peaks agree very well with the predicted widths of envelope solitons calculated 
using the observed values W L /n e T e and other solar wind parameters. These observations 
support the view that strong turbulence processes, namely the modulational instability 
or OTSI, are the means by which some of the type III electron beams are stabilized. 

In some cases, low frequency electric fields show weak spectral enhancements at 
~ 100 Hz. perhaps corresponding to long wavelength ion acoustic waves. The close 
association of these waves with the modulationally unstable Langmuir solitons, indicate 
that the electrostatic decay instability coexists with the modulational instability in 
certain type III burst source regions. 
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Figure 1. Time profiles of the electric field components corresponding to the harmonic and 
fundamental local type IK solar radio bursts and the corresponding Langmuir waves. The 
columns I, 2 and 3 correspond to 11 December 1990, 22 February 1991 and 7 March 1991, 
respectively. The first and second rows correspond to type III bursts at second harmonic 

and the fundamental of the electron plasma frequency, f pe , and the third row corresponds to 
Langmuir waves. 

Figure 2. High time resolution observations of Langmuir waves observed during the type III 
burst events of (a) 11 December 1990, (b) 22 February 1991, and (c) 7 March 1991. These 
measurements were taken with a wide band (6-60 kHz) filter, and 1.12 msec time resolution. 
The enhanced background level after about 600 ms in many of the panels was caused by the 
switch-on of the 32 dB attenuator. Narrow intense spikes observed on top of the central broad 
peaks of the Dec. 11 and March 7 events have been removed from the data. The broad central 
peaks seen on all three days satisfy the conditions expected of envelope solitons. 

Figure 3. Observed 0.2-power width L 0 . 2 of the broad FES peaks versus (n e T e )/W L ) 1 / 2 . Here 
the Lo .2 is calculated using the relation. L 0 . 2 = t 0 2 x V 3W , where t 0 . 2 is the 0.2-power duration of 
the broad peaks and V, w is the solar wind speed as given in Table 1, n e and T e are the electron 
densitv and temperatures, respectively and W L is the peak energy density of the broad FES 
peaks. The correlation coefficient in this case is 0.9. 

Figure 4. Spectral plots of electric fields observed by the YVFA and PFR during 1- or 2-minute 
intervals covering the high time resolution snapshots of Fig. 2. The dotted line corresponds to 
instrumental threshold. The spectral peak in these figures at ~ 10 4 Hz corresponds to Langmuir 


waves. 



Table 1 . Summary of the observed solar wind parameters during the type HI 
events of the present study. Col. 1 gives the date of type III event; col. 2, the ob- 
ed electron density, n e in m 3 , col. 3, the observed electron plasma frequency. 
f P ' in kHz, col. 4, the observed magnetic field, B in units of nT, col. 5, (Q e /uj pe ) 2 , 
where fi e and w pe are the electron cyclotron and electron plasma frequencies, re- 


spectively, col. 6, the average electron temperature, T e , col. 7, the average ion 
temperature, T„ col. 8, the Debye length X D in meters, col. 9 gives the solar wind 
speed, V sw in km/s, and col. 10 gives the sound speed, c, in km/s. 
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Table 2. Summary of the high time resolution observations of Langmuir waves. 
Here E L is the observed peak electric field corresponding to the broad central 
peaks in units of V/m and W L /(n t T e ) are the corresponding normalized energy 
densities. The t 0 2 is the measured 0.2-power duration of the broad central peak in 
units of ms and the L 0 . 2 is the corresponding width, calculated using the formula: 
£ 0.2 = *o .2 x V JW , where V sw is the solar wind speed as given in Table 2. The col. 
(3-6) correspond to the respective FES events of the present study. The 5 0 2 is 
the predicted spatial scale of the envelope soliton 
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